Synchronization of flexible sheets 



Gwynn J. Elfring and Eric Lauga** 
Department of Mechanical and Aerospace Engineering, 
University of California San Diego, 9500 Gilman Drive, La Jolla CA 92093-0411, USA. 

(Dated: August 31, 2011) 

When swimming in close proximity, some microorganisms such as spermatozoa synchronize their 
flagella. Previous work on swimming sheets showed that such synchronization requires a geometrical 
asymmetry in the flagellar waveforms. Here we inquire about a physical mechanism responsible for 
such symmetry-breaking in nature. Using a two-dimensional model, we demonstrate that flexible 
sheets with symmetric internal forcing, deform when interacting with each other via a thin fluid 
layer in such a way as to systematically break the overall waveform symmetry, thereby always 
evolving to an in-phase conformation where energy dissipation is minimized. This dynamics is 
shown to be mathematically equivalent to that obtained for prescribed waveforms in viscoelastic 
fluids, emphasizing the crucial role of elasticity in symmetry-breaking and synchronization. 
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I. INTRODUCTION 



Motile microorganisms swim in a fluid regime where inertia is unimportant and viscous stresses dominate. In this 
limit the flow field due to a swimmer affects the motility of nearby cells [1], a fact which is biologically important 
as microorganisms such as spermatozoa are often found in high-density suspensions [2] . A particular consequence of 
'"^ , these fiuid-based interactions is the synchronization of the flagella of some spermatozoa observed to occur when these 
cells are swimming in close proximity [3-5]. 

G.I. Taylor provided the first quantitative analysis of this phenomenon in his landmark work on the mechanics 
of swimming microorganisms [6]. Using a two-dimensional model of sheets passing sinusoidal waves of transverse 
■ displacement, he showed that the energy dissipated by the fluid due to such motions, for two swimmers a fixed 
distance apart, is minimized if they are in-phase [6]. Subsequent computational studies [7-9] have shown that such 



^ ' synchronization can occur due to fluid forces alone. 



In recent theoretical analysis it was demonstrated that two infinite sheets passing waves of a prescribed shape, 
will not synchronize in a Newtonian fiuid if the shape of the waveforms r^i and •q2 satisfy 772(2;) = —rji{—x -\- 9) 
(where is a fixed phase shift and x is the direction along the sheets) because of the kinematic reversibility of the 
Stokes equations [10, 11]. The sinusoidal waveforms of Taylor's swimming sheet clearly fall into this category, and 
" ' \ will thus not dynamically synchronize in a Newtonian fiuid. It has been observed that excess symmetry similarly 
^ . curbs synchronization in other models [12-16]. For a sinusoidal sheet, a geometric perturbation must therefore be 
' added (for example in the form of a higher order mode) to break the necessary front/back symmetry, and give rise 
f — , to a time-evolution of phase toward the synchronized state [10, 11]. Alternatively, instead of a geometric symmetry 
' breaking, it has also been shown that synchronization can occur if the kinematic reversibility of the field equations is 
0^ \ removed, as is the case for a viscoelastic fiuid [17]. In such a scenario the phase always evolves to a stable in-phase 
■ conformation where the energy dissipated by the swimmers is minimized. 

\ In this paper we inquire about a physical mechanism responsible for symmetry-breaking in real biological cells. 
T-H , Instead of delineating a fixed waveform for the swimming sheets, we take the more realistic modeling approach of 
passing internal waves of bending as produced by a flagellum's internal structure (or axoneme) [18]. In the case of a 
single sheet, similar models have been employed to study swimming [19, 20] and peeling [21]. We use this model to 
' show that elastic deformation due to fiuid body interactions, with purely sinusoidal forcing, always leads to in-phase 
synchronization. 

Flexibility has long been considered as an avenue for symmetry breaking in Stokes fiow. Purcell, in his celebrated 
talk and paper [22] , asserted that while a stiff oar undergoing reciprocal motion would produce no net motion, due to 
the scallop theorem, a flexible oar would escape this conundrum because of a broken symmetry between forward and 
reverse strokes. This was first investigated analytically by Machin [23] while more recently a number of theoretical 
and experimental studies have further elucidated the effect of flexibility quantitively, both for boundary-driven [24-26] 
and internally-driven filaments [27-29]. 

Recently, the analysis of a pair rotating helices as model for bacterial flagella has shown flexibility to be a crucial 
ingredient for synchronization [30], as was similarly shown for a minimal model of interacting cilia [31]. Additionally 
flexibility has been found to be requisite for the synchronization of paddles that would otherwise be too symmetric to 
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FIG. 1. System of two infinite two-dimensional sheets, of shape described by the functions rji and 772, which are separated by 
a fluid layer of thickness h (mean value, h), have a relative phase (j), and may move relative to each other with a velocity Ua- 
An infinitesimal material element on the top sheet is subject to fiuid pressure p, normal force A^, tension T and moment M. 
Inset: schematic representation of in-phase and opposite-phase configurations. 



yield stable fixed points [32]. In these models and experiments, the bodies are rigid, but permitted to deviate from 
their trajectories, in an elastic manner in response to fluid forces. In contrast, here we allow the bodies themselves to 
deform due to fluid stresses induced by the other swimmer. 

Our approach is organized as follows. For a pair of two-dimensional sheets in the lubrication limit, we derive 
a system of nonlinear equations governing both the fluid stresses and the resulting swimmer shapes. We linearize 
these equations to produce analytical solutions, and then solve the nonlinear equations numerically. We show that 
flexible sheets with symmetric sinusoidal forcing will deform when interacting with each other via a thin fluid layer 
in such a way as to break geometrical symmetry, and to evolve to an in-phase conformation where energy dissipation 
is minimized. Further, this evolution of phase is shown to be functionally equivalent to that found for prescribed 
waveforms in viscoelastic fluids [17], illuminating the role of elasticity in symmetry-breaking and synchronization. 



II. MODEL SYSTEM 



We consider the dynamics of two infinite two-dimensional elastic sheets which are separated by a fluid layer of 
mean distance h (see figure 1). The sheets deform due to a balance between an active moment m, passive bending 
(elastic) resistance, and fluid stresses. The positions of the sheets are given by yi = rii{x,t) for the bottom sheet and 
y2 = h + r]2{x, t) for the top sheet. We look to solve this problem in the limit that the fluid layer is thin compared to 
the wavelength of the sheets, kh <^ 1, to make use of Reynolds' lubrication approximation for the fluid field equations 

An infinitesimal element along one sheet is subject to tension, T, normal force, TV and moment, M (shown in 
figure 1). Since the fluid layer is thin and lubrication forces are singular with the gap thickness, forces from the outer 
flow are safely ignored [20] . To capture the active bending of the sheet we use a model for the flagella of eukaryotes 
introduced by Jiilicher and co-workers [18, 28] where the bending of an elastic filament (flagellum) is caused by the 
constrained sliding of microtubule doublets; the effect of the internal forces which induce sliding are here represented 
by an active moment density to. Given that we are in an over-damped limit we take the equilibrium shapes to arise 
instantaneously [21]. Assuming reasonably small deflections such that the swimmers are linearly elastic [34], force (F) 
and moment (M) balance on the top sheet yield respectively 

dF dM 
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where n is the unit normal and <t is the fluid stress tensor. The relation between the moment and the sheet curvature 
is given by the constitutive relation M = Bd^rj/dx^, where B is the sheet bending stiffness. Combining with (1) we 
obtain the equations governing the shape of the sheets 

9^77 dm , , 

where =F are for the top and bottom sheets respectively. If resistive force theory is used for the fluid forces then we 
obtain the governing equation used by [18] for a single filament. 

The internal forcing on the sheets is assumed to take the form 771(2;, €) = Ag{kx — kct) where k is the wavenumber, 
c is the wave speed, A is the amplitude of the moment and g is an arbitrary but 27r-periodic function. Because 
the forcing is periodic, we will assume the shape 77 to also be periodic. We nondimensionalize vertical distances by 
y* = y/h and horizontal distances x* = kx (* indicates a dimensionless quantity). Nondimensionalizing the continuity 
equation we find that if the horizontal velocity is given by u = cu* then the vertical velocity must he v = ecv* where 
t ~ kh. The Stokes equations then yield the lubrication equations to leading order in e 

dp* ^ (Pv*_ ^ = (3) 
dx* dy*^ ' dy* 

du* ^dv* _ ^ 

dx* dy* 

where p* = e^p/^Lo. Forces (per unit depth) are nondimensionalized as /* = fe/fic, while energy dissipation rate per 
unit depth is E* ~ e^E/ pooch. 

In the lubrication limit, e ^ 1, the normal force due to the fluid on the beam is to leading order merely the pressure, 
— n ■ ex ■ n = p. Since the field equations for the fluid yield the pressure gradient, we differentiate (2), and recasting 
the equation in dimensionless form we obtain 



B*^ = A*-^ + ^ (^) 
dx*^ dx*^ dx* ' ^ ' 



where B* = Be^k^/fiut is the dimensionless bending stiffness and A* — Ae^k'^/fiuj is the dimensionless amplitude of 
the active bending moment. This equation allows us to solve for the shape of the sheets, 771^2, and is coupled to the 
fluid field equations through the pressure gradient. We now drop the * for convenience. 

III. ANALYSIS 

Given the form of the forcing we expect post-transient solutions which are functions of a wave variable z = x — t and 
thus we write 77 = 77(2). The top sheet may move relative to the bottom sheet with a horizontal velocity u = Ua, hence 
the boundary conditions for the fluid equations (4), in a frame moving with waveform are given by u(x,yi) = — 1, 
v{x,yi) = — u{x,y2) ~ Ua — 1 and v{x,y2) = —772- Given the above boundary conditions the solution for the 
velocity field is found to be 

u{x,y)==^^iy-yi)iy-y2) + UA- — ^ - 1. (6) 
2dx 7/2 - yi 

If one integrates the continuity equation one finds 

di ""^y^^^^' 

If Ua = then the fiow rate between the sheets is constant. Integrating (7) and exploiting the periodicity of the 
pressure [10], we obtain the equation for the pressure gradient 

dp ^ 6Ua - 12 _ 12^A2/2 _ {6Ua - 12)/2 - UUaJs 

dx " hh^ ' ^ ' 

where the distance between the two sheets is given /)- = 1 + 772 — 771 and Ij ~ J^^ h~^dx and J3 = J^^ y2h^^dx. Then 
the force on the top sheet is given by 

Jo V 2 da; 1 + 772 - 771 / 



. . dp du 



v=y2 
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A. Linear regime: Statics 



We first look at the case where we enforce Ua = (i.e. we fix the top sheet with respect to the bottom sheet) in 
order to determine under which condition a nonzero synchronization force will arise. 

If we assume the dimensionless amplitude of the forcing A to be small, and the dimensionless bending stiffness B to 
be large, then the shape amplitude (or maximum value of the shape) ||ry||oo is expected to also be small. Linearizing 
the pressure gradient for small 77, with = gives 

J«-12(,72-m), (10) 

where we have invoked an integrated conservation of mass, {r]2 ^ "Hi) =0 (angle brackets (), denote the average over 
a period). Our goal is now to determine whether two sheets which are equally and symmetrically forced but with a 
phase shift will synchronize in time. With this in mind we let 92 = 3i(-z + 4>) with gi{z) = cos(z), and we set equal 
for both sheets the bending stiffness B and the forcing amplitude A. The linearized governing equations are then 
given by 

- 12(7^2 - ryi) ^ -v4cos(z), (11a) 

+ 12(7^2 - ?7i) = -^cos(z + ^). (lib) 

To solve these equations we apply periodic boundary conditions. The solution of this system of equations with the 
linearized pressure gradient (10), can be found analytically to be 

12B [ cos(2 + ^) - cos 2] - [ (288 + sin z + 288 sin(z + (/))] 

'^^(^^ = ^ BWmTW) + (^'"^ 

12B [cos z - cos{z + (f)] - [288 sin z + (288 + B'^) sin(z + 0)1 , , 

= ^ BWf^TW) + ^- (^2^) 

Note that nothing prevents the solution from including a uniform shift C{A,B); however, the relevant physics of the 
problem are invariant under such a shift hence and hence the value of C is irrelevant (equivalently, we place our z-axis 
at (7?i) =0). 

Both shapes in (12) are delineated by the competition between bending rigidity, the pressure gradient in the fluid, 
and the internal forcing (with (j) dependence). In this linear limit we note that the sheets are linear in the forcing 
amplitude A, and when A = Q then as expected they become straight, i.e. t] = 0. In the limit where rigidity dominates, 
B — >■ 00, then the sheets also tend to become straight, 77 0. If the rigidity and forcing amplitude are both very large 
{A, i? ^ 1) then we can scale out the contribution from the fluid forces in (11) and we are left with rji « ~{A/B) s'mz 
and 772 ~ sm{z + (f) as might be expected. Note finally that the solutions in (12) are only valid when ^ A 

as otherwise unphysical solutions may arise with the sheets overlapping; this is prevented when the full nonlinear 
form of the pressure gradient is kept, as it diverges when h ^ Q. 

Solutions to (12) for A = 1, i? = 10, and three values of the phase difference {(p — 7r/4, 7r/2, 37r/4) arc plotted in 
figure 2 (left). We observe that the shapes are sinusoidal and hence individually remain symmetric both about the 
vertical axis and the horizontal axis. The global asymmetry that arises however is that the amplitudes of the two 
waveforms are not equally modulated by the fluid pressure, as shown in figure 2 (right). We see indeed that the top 
sheet has smaller amplitude for positive (j) (and by symmetry, this is reversed upon changing (p ^ —(p). 

The phase locking force on the top sheet is, at leading order, given by 

fx = -Q (772 - 771) (?72 + ?7i)d.T = 27ra sine/). (13) 



^0 

where a — 14AA^ / {576B + B^) . Equation (13) is the main result of our paper. The phase locking force is proportional 
to the sine of the phase, meaning that the only stable fixed point is expected to occur at (j) = 0, and hence all initial 
conformations will evolve to the stable in-phase conformation. We thus see that the elasticity of the swimmers, and 
thus fluid-structure interactions, can introduce the geometrical symmetry-breaking necessary to develop a nonzero 
phase locking force. The force is found to be quadratic in amplitude, reminiscent of viscoelastic symmetry-breaking 
[17]; by comparison, the phase-locking force arises at fourth order in amplitude for prescribed asymmetric waveforms 
in a Newtonian fluid. The reason for the difference is that with elastic deformation, the sheets are ultimately not the 
same shape, despite having identical mechanical properties, and hence (77J) ^ {^2)1 prescribed waveforms this is 
different as the same waveform is prescribed for both sheets, and the quadratic term of the force vanishes. 
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FIG. 2. Left: Solution shapes, (12), for various phase differences = 7r/4 (sohd lines), 7r/2 (dashed lines), and 37r/4 (dash-dot 
lines), with A = 1 and B — 10. We observe that the amplitude is not evenly affected by the pressure (the plot is shown here 
with h — 0.2 rather than h = 1 for display purposes only). Right: Shape amplitude, ||?7||oo, vs. phase difference, (bottom 
sheet: solid line; top sheet: dashed line). Lines are reversed upon the change —(p. 



The energy dissipated by the fluid between the two svifimming cells to leading order is 

E = l2j^ im-m) dx= ^^^_^^^ il~cos^). (14) 

We see that the energy dissipation is a global minimum when cj) = and global maximum when (f) = tt. It follows then 
that the cells will always evolve to a state of minimum energy dissipation. We observe that the form of the energy 
dissipation is precisely the same as that for fixed shapes (and taking the fixed wave amplitude A'j^^^^ ~ / {576 + B^) 
they are equal) [11]. It is important to note that waveforms with a prescribed broken symmetry may evolve to either 
the in-phase or opposite-phase conformation [10]; in contrast, the natural symmetry-breaking due to elasticity of the 
bodies, or in the fluid, leads to a conformation of minimum energy dissipation. 



B. Linear regime: Dynamics 



When the sheets are permitted to evolve in time in force-free swimming, the relative velocity Ua will thus be 
nonzero. In order to determine the leading order component of the pressure field we must first find out how the 
relative velocity scales with the sheet amplitudes. Given that the net force on the sheets in the dynamic case must 
now be zero, we obtain at leading order the relative speed as given by 



3 '•2^ 



Ua = / iv2 - m)iV2 + Vi)dx. (15) 

Jo 

We see that the velocity is quadratic in amplitude and indeed is proportional to the static force Ua = /J/Stt (we use 
here the superscript s to indicate the static force given by (13) to avoid confusion). With this we thus know that at 
leading order the pressure field (for a given 0) is invariant between the static and dynamic case, and since the beam 
equation couples via the pressure field, our instantaneous shapes are found to be the same. The only difference is that 
now the phase difference changes in time (geometrically) due to the presence of a nonzero relative velocity, according 
to d(j)/dt = —Ua- Using (15), we find the rate of change of the phase at leading order to be given by 

^ = — asint/). (16) 

Equation (16) can be integrated analytically, leading to a time-evolution of the phase given by 



(j){t) = 2tan" 



tan I ^ 



(17) 



All initial conformations, (j}Q, decay in time to the stable in-phase conformation, = 0. Notably, the time-evolution of 
the phase for a sinusoidally forced elastic sheet we obtain here is mathematically similar to that for a fixed sinusoidal 
waveform in a viscoelastic fluid [17] and for rigid bodies with flexible trajectories [31], emphasizing therefore the 
crucial role of elasticity in synchronization. 
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FIG. 3. Top: solution shapes, rji and ri2, for bending stiffness B = 10, phase difference (p = 37r/4, and amplitudes A = {1, 5, 10} 
(left: Ua = 0; right fx = 0); linearized pressure approximation: solid lines; full pressure gradient: dashed lines. Bottom 
left: phase locking force, fx/A^, vs. phase difference, 0, when Ua ~ 0; bottom right: relative velocity, Ua/A^, vs. (f>, when 
fx ~ 0. Both plots are for numerical solutions of the nonlinear equations with B = 10 and A = 1 (circles), 5 (squares), and 
10 (diamonds); linearized solutions are shown as solid lines. Away from the linear regime, the rate of change of the phase is 
affected while the forces are not. 

Finally, if we allow a small difference in the wavespeeds of the sheets. Aw, then to leading order we have the same 
evolution of phase, (16) and (17), but now the rate of change of phase is defined as (/) = —Ua — Acj and hence we see 
a synchronization of shape but not of material points. 



Nonlinear Case 



To move beyond the linear regime, we now solve the nonlinear equations for the shapes (5) numerically, together 
with (8), using Matlab's boundary value problem solver bvp4c, both for the static case, Ua = 0, and the force free 
case, fx = using (9). 

We find the linearized pressure gradient to be a capable approximation, particularly when the bending is of the 
same order as the pressure _B ~ 1 and the phase difference is small; however, when A, B ^ I the linearized pressure 
may lead to unphysical solutions particularly if the sheets are near opposite-phase as the divergent nature of the full 
form of the pressure gradient is required to deform the sheets from contact. 

In figure 3 we illustrate the breakdown of the linear regime. We plot the static shapes {Ua — 0, top left), and 
dynamic shapes (fx — 0, top right), both with a phase difference oi 4> — 37r/4, bending stiffness B ~ 10, and with 
forcing amplitudes A = {1, 5, 10}. We see that for increasing amplitude the shapes predicted by linearized equations 
(solid) and nonlinear equations (dashed) begin to diverge. In particular the nonlinear equations lead to a pronounced 
left-right asymmetry in the individual shape. In the lower left plot of figure 3 we display the phase locking force 
vs. phase, while in the lower right plot we show the relative velocity vs. phase, both with i? = 10 and A = {1, 5, 10} 
(circles, squares and diamonds respectively) for the numerical solutions to the nonlinear equations; the analytical 
solutions for the linear equations are shown solid. The synchronizing hydrodynamic force, fx, found by either method 
remains remarkably consistent even for very large forcing amplitude, A. In contrast, the rate of change of the phase 
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FIG. 4. Time- evolution of the phase diflerence, in the nonlinear problem for B = IQ starting from an initial angle of 
<j>/n — .99 and with forcing amplitude A = 1 (circles), 5 (squares) and 10 (diamonds); the linear estimate is shown solid. As 
the forcing amplitude increases, the linear solution increasingly underestimates the time scale to synchronize. 

decreases markedly from the linear approximation for large amplitude waves. Because the force is virtually unaffected 
we know therefore that the resistance to motion is dramatically increased by the change in shape. 

In figure 4 we integrate the instantaneous relative velocity to obtain the evolution of the phase in time. We show 
solution to both the linear equations (solid) and nonlinear equations (dashed) for B = 10 and A = {1, 5, 10} starting 
from an initial phase difference of just less than tt. We see that, as the forcing amplitude increases, the nonlinear 
equations yield an increasingly slower evolution to a synchronized conformation than that predicted by the linear 
regime; however, the general behavior remains qualitatively unchanged. 

IV. CONCLUSION 

In this paper we inquired about a physical mechanism responsible for symmetry-breaking and synchronization in 
the flagella of biological cells such as spermatozoa. In a Newtonian fluid, two swimming sheets passing waveforms of 
a prescribed sinusoidal shape will not synchronize due to an excess of symmetry; however, here we have demonstrated 
that identical flexible sheets with symmetric sinusoidal forcing will deform, when interacting with each other via a thin 
fluid layer, in such a way as to systematically break the overall geometrical symmetry. This system will always evolve 
to an in-phase conformation where energy dissipation is minimized, in contrast to a prescribed asymmetry, which 
may maximize energy dissipation. In addition, this time-evolution of the relative phase is shown to be equivalent 
to that obtained for prescribed waveforms in viscoclastic fluids [17], emphasizing the crucial role of elasticity in 
symmetry-breaking and synchronization - be it that of the fluid, or the swimmers themselves. 
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